function [eps] = randeps(N,amp)
% randeps(N,amp)
% 	Parameters:
% 		N:	number of points in the time series
% 		amp:	the maximum value of abs(eps) for each epsilon to be generated
% 	Outputs:
% 		eps	a sequence of values of length N where each value is between -amp and +amp.  
% 	Purpose of Code:  This code is designed to produce a sequence of values length N where each value is randomly selected from a set of values
% 	that are uniformly distributed over the interval [-amp, amp].  This code is a helper for the code twoshiftpddyneps.

for j = 1:N
    eps(j) = ((2.*rand(1))-1).*amp;
end


end
